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The  power  spectra  of  P  wave  signals  are  studied.  Data  from  the  ?2  long 
period  (LP)  and  short  period  (SP)  sensors  of  the  NORSAR  array  in  Norway 
xs  used.  For  each  type  of  sensor,  the  signals  from  all  sensors  are  com¬ 
bined  in  an  optimum  way  to  produce  average  LP  and  SP  spectra.  Bias  due 
to  noise  power  is  removed,  the  instrument  responses  are  corrected  for 
and  the  LP  and  SP  spectra  merged  together.  About  150  events  are  treated 
of  these  the  LP  signals  are  analyzed  for  about  20  events.  The  earth¬ 
quake  and  explosion  spectra  agree  well  with  the  Dredictions  of  analvti- 
cal  source  models.  The  corner  frequencies,  where  the  power  begins  to  de 
crease  with  increasing  frequency,  are  roughly  three  times  higher  for 
explosions  than  for  earthquakes  of  comparable  magnitude  (in  the  mp 
range  5. 0-6.0).  This  causes  the  low-frequency  power  level  to  be  lower 
for  explosions.  At  the  frequency  0.3  Hz  explosion  siqnal  power  is  typi¬ 
cally  20  dB  less  than  the  power  of  signals  from  earthquakes  with  simi¬ 
lar  magnitude.  Large  and  systematic  regional  variation  is  observed  in 
the  slope  of  spectra  at  frequencies  above  2  Hz.  This  is  interpreted  to 
stem  from  variation  of  anelastic  attenuation  in  the  Earth.  This  large 

VJr:ut;Lu?  uake!  suspect  the  value  of  event  classifiers  employing  energy 
a  the  higher  frequencies  of  the  P  wave  spectrum.  Comparisons  of  power 
in  the  vicinity  of  the  corner  frequency  may  offer  classification  capa- 
oxlity  because  the  expected  spectral  differences  are  well  demonstrated 
and  the  influence  of  attenuation  during  wave  transmission  is  not  yet 
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FOREWORD 


The  work  reported  in  this  document  was  performed  at  Insti¬ 
tute  of  Seismology,  University  of  Helsinki  under  the  grant 
OSR-72-2377 .  The  author  visited  the  Norwegian  Seismic  Array 
Data  Center  in  Kjeller,  Norway  during  the  first  half  of  the 
period  covered  by  this  report,  and  the  significant  help  re¬ 
ceived  from  the  NORSAR  staff  is  gratefully  acknowledged. 

The  author  could  benefit  from  innovative  discussions,  in  par¬ 
ticular  with  the  staffs  of  the  Seismic  Array  Analysis  Center 
ii:  Alexandria,  Virginia  and  of  the  Lincoln  Laboratory,  MIT, 
Boston,  during  a  short  trip  to  the  United  States.  Thanks  are 
due  to  several  very  helpful  persons  met  during  that  visit. 
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I .  INTRODUCTION 


In  this  work  the  power  spectra  of  P  waves  from  seismic  events 
recorded  at  the  NORSAR  large  seismic  array  in  Norway  have  been 
observed.  An  objective  of  the  study  is  to  compare  the  shapes 
of  spectra  from  earthquakes  and  explosions,  in  the  hope  of 
finding  stable  differences  between  the  two  types  of  events, 
which  would  aid  in  the  classification  of  seismic  events.  Other 
objectives  are  to  estimate  parameters  of  the  seismic  sources 
by  comparing  the  observed  spectra  to  predictions  made  using 
simplified  source  models1'2'3'4'5'6  and  to  observe  regional 

variations  of  spectra  which  could  be  attributed  to  propagation 
path  effects. 


In  this  report  observations  of  P-wave  spectrum  in  the  frequen¬ 
cy  range  0.03  and  6  Hz  are  described.  Observations  in  a  broad 
frequency  range  are  also  reported  in  references7,8.  Wyss,  Hanks 
and  Liebermann  reported  the  greatest  differences  between  earth¬ 
quake  and  explosion  spectra  to  be  found  at  low  frequencies, 
where  the  signal-to-noise  ratio  (SNR)  for  weak  events  is  ooor . 

For  classification,  differences  in  the  frequency  range  0.5-4  Hz 
which  has  a  good  SNR  would  be  most  valuable.  That  frequency 
range  has  been  used  for  computing  single  numbers ( ’discriminants ' 
which  describe  the  relative  levels  of  the  high  and  low  frequency 
energy  of  a  spectrum9'10'!1' 12 ■ 13' I4 .  In  thls  VJr*  such  a  sin_ 

ql°  Parameter  is  also  defined  and  computed  for  a  number  of  spec¬ 
tra.  The  value  of  such  parameters  for  classification  is  still 
unclear. 

In  estimating  the  spectrum  of  a  P  wave  using  a  large  seismic 
array  an  optimum  method  for  combining  data  from  all  sensors 
should  be  used.  The  NORSAR  short  period  sensors  give  132  paral- 
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lei  recordings  from  each  signal.  Each  of  these  differs  slightly 
from  all  the  others,  while  the  noise  at  the  different  sensors 
is  assumed  to  be  completely  uncorrelated.  To  make  use  of  all 
sensors  and  to  arrive  into  a  single  power  spectrum,  a  basic  de¬ 
cision  is  whether  one  should  average  over  the  signals  them¬ 
selves  {beamforming)  prior  to  computing  the  power  spectrum,  or 
whether  the  average  over  the  power  spectra  computed  for  each 
single  sensor  should  be  used  (spectraforming) . 

Lacoss  and  Kuster15  showed  that  (a)  the  spectra  computed  by 
spectraforming  and  beamforming  differ  from  each  other  ever  in 
the  absence  of  noise  and  that  (b)  in  case  of  weak  signals  the 
selection  of  a  proper  method  is  of  importance  in  extracting  the 
signal  from  the  noise.  The  delays  required  for  the  delay-and-sum 


i  ’  ’I’  .N0RfAR  array  configuration .  The  22  subarrays  are 

shown.  Each  subrray  is  equipped  with  6  vertical  short-period 
sensors  and  one  3-component  long-period  seismometer.  ? 
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operation  in  beamforming  may  be  poorly  known  for  a  weak  event 
and  this  may  cause  unpredictable  and  large  beamforming  losses 
at  the  higner  signal  frequencies.  This  is  a  strong  argument 
against  the  use  of  beamforming  in  computing  power  spectrum  at 
the  the  higher  frequencies.  If  no  beamforming  losses  would  ex¬ 
ist,  beamforming  would  be  the  best  method  for  extracting  the 
signal  spectrum  from  the  noise,  but  tne  losses  worsen  its  per¬ 
formance  . 

All  spectra  presented  in  this  report  are  computed  using  a  com¬ 
bination  of  beamforming  and  spectraforming.  We  averaged  22 
power  spectra  computed  for  the  22  subarray  beams  from  the  NORSAR 
subarrays.  Each  subarray  beam  was  formed  by  averaging  signals 
with  suitable  delays,  recorded  by  6  sensors  within  a  region  of 
7-8  km  in  diameter.  The  geometry  of  NORSAR  is  shown  in  fig.l.l. 

In  terms  of  stability,  beamforming  loss  and  computing  time 
this  was  considered  a  good  compromise.  In  Appendix  1  measure¬ 
ments  of  SNR  and  beam  losses  from  events  in  Central  Asia  are 
presented,  and  the  reasons  for  the  selection  of  the  procedure 
of  averaging  subarray  beams  are  presented  in  more  detail.  These 
measurements  have  been  published  earlier16. 

In  extracting  the  spectraformed  signal  from  the  noise  the  most 
important  operation  is  the  subtraction  of  an  estimate  of  noise 
power  from  the  signal  power.  With  certain  assumptions  on  noise, 
the  expected  value  of  the  remainder  is  equal  to  true  signal  pow¬ 
er  in  the^ frequency  band  considered15.  As  noted  also  by  Linde 
and  Sacks  ,  this  removal  of  bias  due  to  noise  is  of  prime  im¬ 
portance  when  observing  frequencies  out  of  the  dominant  signal 

frequency.  The  noise-correction  is  described  in  detail  in  Appen¬ 
dix  2 . 

In  collecting  data  it  was  tried  to  include  events  from  various 


tectonic  provinces  and  to  observe  earthquakes  in  the  vicinities 
of  nuclear  test  sites.  All  included  earthquakes  were  shallow, 
as  a  rule  maximum  depth  limit  of  50  km  was  used,  but  including 
some  events  in  the  depth  range  50  to  70  km.  The  depth  70  km 
corresponds  to  the  supposed  lower  boundary  of  the  low-attenua¬ 
tion,  high-Q  lithosphere.  Initially  only  events  in  the  magni¬ 
tude  range  5.0  -  6.0  m^  were  analyzed,  because  a  homogeneous 
data  set  was  desired.  That  also  is  the  optimum  observation 
range  if  both  short  and  long  period  sensors  are  used.  Later  a 
wide  magnitude  range  was  covered.  A  list  of  all  included  events 
is  given  as  Appendix  3. 

The  magnitudes  of  events  quoted  in  this  report  are  P  wave  mag¬ 
nitudes  (m^)  given  by  the  (U.S. ) National  Earthquake  Information 
Center  (NEIC) ,  formerly  referred  to  with  abbreviations  USCGS 
and  NOAA.  Whenever  the  magnitude  is  given  from  another  agency 
or  is  computed  from  other  waves,  this  is  indicated. 

In  computing  a  power  spectrum  using  long  period  sensors,  signals 
from  the  22  vertical  LP  sensors,  one  in  each  subarray,  are  beam- 
formed  (summed  with  delays) ,  a  signal  section  96  seconds  long 
and  beginning  1-2  sec  prior  to  the  onset  time  of  the  P  waves  is 
Fourier  transformed  and  the  absolute  value  of  the  transform  is 
squared.  Bias  due  to  noise  is  removed  by  subtraction  of  noise 
power.  The  spertra  presented  in  the  figures  of  this  report  are 
smoothed  over  a  bandwidth  of  0.021  Hz.  Tn  computing  a  power 
spectrum  using  3hort  period  sensors,  the  combined  bcamforming/ 
spectra forming  procedure  described  above  was  used.  Signal  sec¬ 
tions  of  6.4  or  12.8  seconds  were  used.  The  noise-corrected 
power  spectra  presented  in  figures  were  smoothed  over  a  band¬ 
width  of  0.313  Hz. 
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II.  GENERAL  SHAPE  OF  SPECTRA 


In  this  chapter  the  general  shape  of  earthquake  and  explosion 
spectra  is  described  using  spectra  computed  from  shallow  earth¬ 
quakes  and  explosions  in  Eastern  Europe,  Western  and  Central 
Asia,  at  the  western  coast  of  North  America  and  at  the  Aleutian 
Islands.  By  combining  measurements  made  with  long  and  short  pe¬ 
riod  (LP  and  SP)  sensors  the  shapes  of  the  P  wave  spectra  with¬ 
in  the  frequency  range  0.03  to  6.0  Hz  could  be  studied.  In  sev¬ 
eral  cases,  however,  the  signal-to-noise  ratio  (SNR)  was  not 
sufficient  below  the  frequency  0.3  Hz  for  the  signal  spectrum 
to  be  extractable  from  the  noise.  In  addition  to  the  observa¬ 
tions  reported  in  this  chapter  for  a  number  of  events  only  the 
signals  recorded  by  the  SP  sensors  were  spectral  analyzed. Those 
observations  will  be  described  in  the  next  chapter. 

In  fig.  II. 1  are  shown  displacement  power  spectra  for  6  earth¬ 
quakes  in  Central  Asia.  The  spectra  are  corrected  for  the  NORSAR 
SP  and  LP  instrument  responses.  The  spectra  shown  are  selected 
from  a  suite  of  13  spectra  as  those  which  had  the  best  SNR  in 
the  long  period  range.  This  will  give  a  somewhat  biased  view, 
but  a  low  signal  power  was  not  the  only  reason  for  rejecting 
spectra.  The  large  variation  in  the  background  noise  level,  in 
some  cases  caused  by  interfering  events,  also  influenced  the 
detection  of  LP  P  waves. 

The  general  impression  is  that  for  these  events  having  magni¬ 
tudes  5.2  -  6.0  (m^)  the  spectra  follow  the  prediction  of  Aki4. 
The  power  level  is  roughly  constant  at  frequencies  below  0.2  - 
1.0  Hz,  at  higher  frequencies  the  power  decreases  rapidly  and 
almost  monotously  with  increasing  frequency.  The  variations  in 
the  relative  levels  of  high  and  low  frequency  energy  are  large. 
The  shocks  101,  102,  105  and  111  all  have  the  magnitude  5.6(n,  )• 
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Fiq.  II. 1,  Earth  displacement  power  spectra  of  earthquakes 
in  Central  Asia.  The  spectra  are  displaced  by  steps  of  15  dB 
for  clarity.  Short  horizontal  lines  at  left  margin  show  for 
each  spectrum  the  level  of  55  dB  over  power  density  1  nm^/Hz 
Magnitudes  are  given  in  parentheses. 


Thi >  is  reflected  in  their  equal  levels  of  power  at  the  frequen¬ 
cy  -.0  Hz,  where  the  measurement  for  the  magnitude  is  made. 
But  their  long  period  power  levels  var  by  20  JI3,  or  by  one  or¬ 
der  of  magnitude. 

This  vai iation  is  coupled  with  the  variation  of  the  corner  fre¬ 
quency,  where  the  spectral  power  begins  to  fall  off  with  in¬ 
creasing  frequency.  It  ranges  from  0.2  to  1.0  Hz  for  these  spec¬ 
tra.  if  the  conclusion  of  Brune5 (1970) ,  that  the  corner  frequen¬ 
cy  of  s  wave  spectrum  is  inversely  proportional  to  the  radius  of 
the  source  fault,  can  be  applied  to  P  waves  the  source  radii  va¬ 
ry  by  a  factor  of  five,  even  within  the  four  shocks  with  magni¬ 
tude  5.6. 


In  fig.  II. 2  the  effect  of  non-elastic  attenuation  is  shown  for 
vinous  values  of  the  frequency-independent  quality  factor  Q. 


The  curves  show  the  magnitude  of  a  factor  exp(-7tf£)  with  t*  as 
the  parameter.  The  quantity  t*  is 


and 


(2.1) 


where  t  is  the  travel  time  and  5  is  the  apparent  average  Q 
along  the  rav  path.  The  available  scanty  evidence  or  t*  along 
the  path  Central  Asia-NORSAR  suggests  a  value  less  than  0.318'19 
Any  value  of  t*  between  0  and  0.3  does  not  change  the  conclu¬ 
sions  on  the  existence  and  location  of  the  corner  frequencies  in 
fig.  II. 1. 


,s  srrSJfr- 

tn°pLKrSlty  1  ""  /HZ-  Ha*"itudes  °f  ^  are^given 
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In  fig.  II. 3  are  presented  spectra  from  three  shallow  earth¬ 
quakes  at  the  western  coast  of  North  America  in  the  magnitude 
range  4.8  -  5.6  n^.  The  spectra  have  a  somewhat  more  peaked 
character  than  the  Central  Asian  spectra.  It  has  to  be  noted, 
however,  that  the  total  response  of  the  earth  to  a  wave  sent 
from  a  source  modifies  the  origina..  spectrum,  particularly  in 
case  of  the  long  period  waves.  An  example  of  this  is  the  inter¬ 
ference  of  the  pP  and  PcP  waves  with  the  P  wave.  In  case  of 
these  events,  both  of  these  later  phases  arrive  within  the 
96  sec  long  window  of  analysis. 

The  corner  frequencies  vary  from  0.2  to  0.5  Hz.  The  slope  of 
spectra  at  high  frequencies  is  clearly  larger  than  that  of  the 
Central  Asian  spectra.  This  is  most  probably  due  to  a  difference 
in  the  total  effect  of  the  non-elastic  attenuation  along  the 
two  ray  paths. 

In  fig.  II. 4  the  spectra  cr  9  explosions  in  Eastern  Kazakh  are 
shown.  No  signal  power  could  be  extracted  from  the  LP  signals 
below  the  frequency  0.3  Hz.  Whenever  observable,  the  power  at 
frequencies  below  1.0  Hz  is  rather  constant.  Comparison  with 
fig.  II. 2  shows  that  if  these  spectra  were  corrected  for  atten¬ 
uation  taking  t  to  be  0.2,  in  all  cases  the  power  would  be 
constant  below  1.0  Hz  or  have  a  maximum  around  that  frequency. 
The  value  t*  =  0.2  corresponds  to  Q  =  2190. 

Seggern  and  Blandford6  predicted  the  far-field  P  wave  spectrum 
from  underground  explosions  to  be  flat  below  a  corner  frequen¬ 
cy  depending  on  the  yield,  with  a  possible  overshoot  at  the  cor¬ 
ner  frequency.  The  picture  emerging  here  is  not  in  disagreement 
with  their  prediction.  The  frequencies  of  the  maxima  observed 
range  from  0.5  to  1.0  Hz  (events  512,  514  and  516)  which  fits 
well  with  the  prediction  of  Seggern  and  Blandford  for  explosion 
yields  50-100  kt.  These  events  have  been  estimated  to  lie  with¬ 
in  the  yield  limits  20  to  200  km20. 


Fig  II. 4.  Earth  displacement  power  spectra  of  underground 
explosions  in  Eastern  Kazakh.  The  spectra  are  displaced  in 
steps  of  15  or  25  dB  for  clarity.  Short  horizontal  lines  at 
left  margin  show  for  each  spectrum  the  level  of  15  dB  above 

power  density  1  nm  /Hz,  Magnitudes  of  events  are  given  in 
parentheses.  y  1 
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One  would  apparently  prefer  to  have  observations  at  lower  fre¬ 
quencies  to  make  certain  the  detection  of  a  corner  frequency. 
There  is  unfortunately  an  inherent  difficulty.  As  soon  as  the 
increase  of  signal  power  with  decreasing  frequency  begins  to 
level  off,  the  signal  is  overwhelmed  by  tne  even  more  steeply 
increasing  noise  and  becomes  impossible  to  measure,  even  with 
the  capabilities  of  a  large  array.  Further,  when  tr.he  yield  in¬ 
creases,  the  corner  frequency  decreases  coming  to  a  region  of 
worse  SNR,  which  offsets  the  general  improvement  in  the  SNR. 

These  difficulties  are  somewhat  overcome  in  observing  the  spec¬ 
tra  of  the  presumed  explosions  in  Western  Russia,  because  they 
lie  closer  to  NORSAR,  at  distances  15°-20°.  In  fig.  II. 5  are 
shown  the  spectra  of  7  events  in  Western  Russia  with  magnitudes 
ranging  from  3.2 (given  by  NORSAR)  to  5.5.  The  larger  events  of 
this  set  have  relatively  more  energy  at  low  frequencies  than 
the  smaller  events.  LASA  array  locates  the  event  WR2  to  Western 
Kazakh,  rather  than  to  Western  Russia.  It  is  thus  not  compara¬ 
ble  with  the  other  events.  Its  LASA  magnitude  is  3.8. 

In  case  of  the  events  WR3  and  WR4  with  magnitudes  4.6  and  4.5 
(as  given  by  LASA)  there  is  a  very  slow  decrease  of  power  with 
increasing  frequency  in  the  observation  range  0.5  to  3.3  Hz. 
Taking  for  these  ray  paths  in  the  upper  mantle  a  Q  less  than 
1500,  a  correction  for  the  effect  of  attenuation  will  result 
into  virtually  constant  power  levels  in  the  observed  frequency 
range.  This  would  mean  the  corner  frequencies  to  be  above  2  Hz. 

The  three  events  with  magnitudes  above  5.0  have  spectra  which 
even  with  a  reasonable  correction  for  attenuation  display  cor¬ 
ner  frequencies  around  or  below  1.0  Hz.  This  is  in  agreement 
with  the  observations  from  E. Kazakh  explosions.  The  event  WR1 
with  NORSAR  magnitude  3.2  shows  a  maximum  of  power  at  2.0  Hz. 


lig.  II. 5.  Earth  displacement  power  spectra  of  presumed  under¬ 
ground  explosions  in  Eastern  Europe.  The  spectra  are  displaced 
in  steps  of  10  dB  for  clarity.  Short  horizontal  lines  at  left 
margin  show  an  equal  power  level  for  each  event.  Magnitudes 
of  events  are  given  in  parentheses. 


lor  the  purpose  of  observing  the  explosion  spectrum  at  frequan- 
cies  well  below  the  corner  frequency,  long  period  P  signals  from 
the  Cannikin  explosion  with  magnitude  6.8  at  Aleutian  Islands 
were  spectral  analyzed  (event  code  122).  In  fig.  II. 6  its  spec- 


strum  is  shown  with  the  spectrum  of  an  shallow  earthquake  at  a 
nearly  same  location.  The  LP  signals  from  these  two  events  are 
shown  in  fig.  n. 7.  The  explosion  is  1-°  magnitude  units  larger 
*"'han  those  Presented  in  previous  figures  and  it  seems  that  the 
corner  frequency  is  between  0.2  and  0.3  Hz. 


Fig.  1 1. 6.  Earth  displacement  power  spectra  for  an  under- 
ground  explosion  at  Aleutians  (code  122)  and  f£  a  la^ge 
earthquake  at  nearly  same  location  (code  123)  The  earthau*k#> 

ITnlsTt  lefthraral„Ushby  ?°  ®  ciarity.  hSSSSl 

tr i  ieft  mar^in  show  for  each  spectrum  the  level  of 

aivAn  ?bOVe  POWwr  density  1  nm  /Hz .  Magnitudes  of  events  are 

cIuHah1?  ^arentheses*  The  lowest  frequency  band  is  not  in¬ 
cluded  in  frequency  smooth! ig.  oc  ln 


This  ranqe  fits  excellently  with  the  corner  frequency  0.26  Hz 
which  has  been  predicted6  for  an  explosion  of  5  Mt  (the  proba 
ble  yield  of  this  explosion20). 
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(code1 122)  andeaLlarqenear^aeakth?  C^nikin  explosion 
The  continuous  lines  represent  h~  Aleu“ans  <code  123). 
events,  the  dots  around^hp^-  C  L?  beams  aimed  towards  the 
recordings.  S  and  N  stand  for  si^ai^d^ 9  Single  cha™el 
preceding  noise,  respectively  Th£  L  2?  f°5  a  samPle  of 
larger  than  the  exDlocHnn  ea,rt^gua}ce  signal  is 

Note  thf  good  coherence  of  i„°ni the  ori9inal  records, 
large  (event  123)?  LP  Signals  w*en  S/N  ratio  is 
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The  decrease  of  power  with  decreasing  frequency  could  be  en¬ 
tirely  caused  by  the  crustal  transfer  function  at  the  source. 

The  transfer  function  is  dominated  by  the  P/pP  interferance. 
Scaling  the  transfer  function  for  a  shallow  explosive  source 
in  a  continental  crust  computed  by  Fuchs  to  the  present 
source  depth  2  km,  we  would  expect  the  power  at  0.05  Hz  to  be 
20  dB  less  than  at  0.35  Hz.  This  is  sufficient  to  explain  the 
slope  of  the  present  spectrum,  and  there  is  no  objection  to 
a  flat  explosion  spectrum  at  low  frequencies  predicted  when 
the  effect  of  the  free  surface  is  absent6. 

The  large  ripples  at  varying  frequencies  superposed  on  the  gen¬ 
eral  trend  of  the  explosion  spectra  in  figs  II. 4  and  II. 5  are 
also  explained  as  successive  bands  of  constructive  and  destruc¬ 
tive  interference  between  P  and  pP.  The  greater  depth  of  the 
earthquakes  causes  these  bands  to  be  narrower  in  the  earthquake 
spectra  and  they  are  smeared  out  in  the  frequency  smoothing. 

The  pP  phase  may  also  remain  outside  the  time  window  taken  for 
analysis. 

In  fig.  II. 8  the  corner  frequencies  read  from  the  earthquake 
and  explosion  spectra  shown  in  figures  II. 1,  II. 4  and  II. 5  are 
plotted  as  function  of  magnitude.  Earthquake  corner  frequencies 
are  roughly  three  times  smaller  than  those  of  explosions.  The 
prediction  for  explosion  corner  frequency  as  a  function  of 
yield  made  by  Seggern  and  Blandford6  is  also  plotted  to  the  fig¬ 
ure,  assuming  the  yield  Y  and  seismic  magnitude  being  con¬ 
nected  by  equation 

log  Y  =  m^  -  3.9 

as  seems  reasonable  for  explosions  in  hard  rock33.  Considering 
the  uncertainties  in  reading  the  corner  frequency  from  a  spec¬ 
trum  and  in  the  magnitude-yield  relationship,  the  fit  is  not 
bad  and  indicates  that  corner  frequencies  are  really  observed. 
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cy  forSexPlosions°in°hardGrock^le 


HI.  SPECTRA  FROM  DIFFERENT  REGIONS 


In  this  chapter  the  shapes  of  spectra  from  events  at  different 
regions  are  compared.  Only  data  recorded  by  short  period  sensors 
is  presented  because  of  the  often  poor  SNR  of  P  waves  at  low 
frequencies.  Because  different  regions  were  observed  to  have 
signifi  'ar.tl],  different  shapes  of  spectra,  the  events  are 
grouped  geographically. 

The  presented  spectra  are  not  displacement  spectra,  as  in  the 
previous  chapter,  but  acceleration  spectra.  Response-corrected 
power  at  each  frequency  has  been  multiplied  by  co4 ,  and  the  pow¬ 
er  scale  in  figures  is  in  dB,  or  equivalent  to  10 *log1Q (P/PQ) 
where  P  is  the  spectral  power  and  Pq  is  some  constant  reference 
level.  The  spectra  within  each  geographical  group  are  normalized 
to  equal  total  power  (at  the  SP  sensor  system  output)  and  the 
spectra  are  averaged  as  described  in  app.  2.  The  individual 
spectra  and  the  average  are  shown  in  figures.  Frequency  smooth¬ 
ing  over  a  bandwidth  of  0.31  Hz  is  done. 

The  events  in  this  chapter  are  usually  identified  only  by  indi¬ 
cating  the  event  set  to  which  they  belong.  Information  on  indi¬ 
vidual  events  arranged  according  to  the  event  sets  is  given  in 
app.  3.  Figures  III.l  and  III. 2  indicate  the  locations  of  the 
events  and  event  sets.  All  earthquakes  included  to  analysis  are 
shallow, i.e.  have  reported  depths  less  than  70  km,  and  most  have 
reported  depths  less  than  50  km.  In  some  cases  the  depth  is  dif¬ 
ficult  to  verify,  but  these  events  are  few  and  generally  at  re¬ 
gions  where  intermediate  or  deep  shocks  are  not  known  to  occur. 

The  15  events  scattered  over  Eastern  Asia,  from  Kamchatka  to 
hin,'> ,  show  a  large  variation  of  spectral  shapes  (fig.  III.  3). 


17 


2l;.III'1'J.An^lyzed  events  excluding  those  in  Europe.  The 
events  are  divided  into  following  event  sets  (numbers  refer 

^  1^Ure  :  EASTERN  ASIA,  2.  SAKHALIN,  3.  KURILES  4  SOUTH 

?F  ??2dUU'  5’  CENTRAL  ASIA,  6.  EASTERN  KAZAKH  (explosions) 

7  GIBRALTAR,  8. WEST  COAST  OF  NORTH  AMERICA,  9.  NEVADA  (ex- 

APT°LA^2)RI  E  U^H/^2TA'  11  •  GULF  °F  CALIFORNIA,'  12.  m»- 
ATLANTIC  ridge.  The  different  symbols  indicate  relative 

gh  frequency  contents  of  spectra,  estimated  by  the  LTM 

measure.  When  several  events  fall  to  a  nearly  same  location 

0  O^TM^nl96  5J™1Valr/S  plotted*  Triangle^  LTM<0 ,  square: 

0 . 0 < LTM ^0.4,  circle:  0.4<LTM<0.7,  asterisk  LTM>0.7. 
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4.  1  ' ,  ,ents  in  Europe.  The  events  are  divided 

lumbers  refer  to  the  figure):  1.  WESTERN  R 
explosion),  2.  WESTERN  KAZAKH  (presumed  explosion)  3 
Sf  CAUCASUS ,  5.  ITALY,  l  YUGOSLAV?^  7 
fere^  symbols  indicate  relative  high-frequency 
spectra,  estimated  by  the  LTM  measure.  Open  circle:  LT] 
circle:  0.4<LTM<0.9,  asterisk  LTM>0.9, 


The  r.m.s.  deviation  of  the  individual  spectra  (in  dB)  from  the 
average  surpasses  8  dB.  On  the  contrary,  when  a  comparable  num¬ 
ber  of  events  is  selected  from  a  small  source  region,  the  spec¬ 
tral  shapes  vary  less.  The  events  in  very  small  source  regions 
at  Sakhalin (fig. Ill . 4 ) ,  Kuriles (fig. III. 5)  and  Japan (fig. Ill . 6) 
each  show  rather  homogeneous  spectra,  the  r.m.s.  deviations 
from  the  respective  averages  generally  being  below  3  dB. 


^ ’u*  Normalized  acceleration  power  spectra  (aster- 

FASTORN  iSTr?ragG  (?frSleS>  °f  thS  15  ShOCkS  °f  the  set 
ef  to  78  degrees?39  5'0'6'3  "b  and  distance  range 


The  spectra  recorded  from  SAKHALIN  differ  clearly  from  the 
KURILES  spectra  in  their  lower  high  frequency  energy  contents. 
One  possible  explanation  for  the  difference  is  a  difference  in 
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e  noneidstic  attenuation  along  the  respective  wave  paths. 
Such  variations  can  explain  why  the  EASTERN  ASIA  set  with  wide 
ly  dispersed  events  displays  a  larger  variation  of  spectral 
shapes  than  the  more  localized  event  sets.  The  magnitude  also 
influences  the  relative  amounts  of  high-  and  low-frequency  en¬ 
ergy,  as  is  observed  for  earthquakes  in  figs  III. 4  and  III. 5 

(SAKHALIN,  KURILES)  and  for  underground  explosions  in  figs 
HI. 6  and  III. 7. 


isks)  withaverages^circlesj^of  ^uakes^n  S**"  <aSter' 

mTf-at^  3?*“  — d  «^?naak1oi!;mtrhLiurf^t46  5 

shocks  in  magnUuden?an?e  m  $!?RS?n'  The  1° 

magnitude  range  5. 0-5. 9  1  V  (2)  The  10  sh°cks  in 


j  i  *  IZ?  ’  Normalized  acceleration  power  spectra  (aster- 

sStSi„p  TVr3965  JCirCleS)  '  (1)  The  10  shocks  of  the  evenl 
set  KURILES  in  magnitude  range  4. 3-4. 9.  (2,  The  2  shocks  of 

the  event  set  KURILES  with  magnitudes  5.7  and  6.0.  The  KURILI 

events  are  all  within  a  50  km  radius  from  44.7  N,  149.3  E, 

at  a  distance  of  69  degrees  from  NORSAR.  (3)  The  9  events  of 

e  event  set  SOUTH  OF  HONSHU  in  the  magnitude  range  3.9-5  2i\ 

The  events  are  within  a  50  km  radius  from  33.2  N? Ul.o  E, 

at  a  distance  of  78  degrees  . 


Normalized  acceleration  power  spectra  (aster- 

sit  EASTERN  KA7?™®.  C  ?!’;  (1)  3  exPl°sicms  of  the  eve. 
r>\  EASTEEN  KAZAKH  in  magnitude  range  3.9(NORSAR)  -  4.4  m,  , 

lalat  r?l0S  °SS  °f  the  Same  event  set  in  the  magnitude  ^ 
’  5*8A;  EASTERN  KAZAKH  is  a  small  source  regior 
f  f  distance  of  38  degrees  from  NORSAR.  (3)  The  15  sh^ks 

SY?6  CENTRAL  ASIA  in  the  magnitude  range  5.0  - 

e.O  m^  and  distance  range  36-56  degrees!  9 


The  explosions  at  different  sites,  even  when  within  similar  mag 
nitude  limits,  show  differences  in  the  relative  amounts  (figs 
III. 6-8  and  III. 12)  of  high-frequency  energy.  These  can  be  ex¬ 
plained  as  differences  in  the  nonelastic  total  attenuation  at 
the  respective  ray  paths.  Explosions  in  Western  Russia (fig. 

III. 7)  display  a  minimum  at  frequency  1.5  Hz,  which  appears  as 


f,1*'7'  Normalized  acceleration  spectra  (asterisks) 

set  WESTE^eRUSSIACl<m  ?f  pre®umfd  explosions  in  the  event 

3 . 2  (NORSAR)  - f  1‘ \  (l)„3„?Vfnts  in  the  ^gnitude  range 
o.zimjKt»ARj  -  4.5  m.  and  distance  ranqe  17-20  e>\  i 

events  in  the  magnitude  range  5  1-5  I  m  (2)  3 

18-25  degrees.  9  3,1  b*3  and  distance  range 


a  minimum  also  in  displacement  spectra.  The  explosion  spectra 
can  be  compared  to  earthquake  spectra  at  certain  regions  where 
explosions  and  earthquakes  occur  at  relatively  clo^e  distances. 
Comparison  between  EASTERN  KAZAKH  explosions  and  CENTRAL  ASIA 
earthquakes  indicates  clearly  differing  spectral  shapes  (fig. 
III. 6).  This  is  disturbed  somewhat  by  the  single  highly  anom¬ 
alous  earthquake  in  Sinkiang  (code  229  in  app.  3).  That  event 
is  not  included  into  averaging  over  earthquake  spectra  in  fig. 
III. 6.  It  will  be  further  discussed  later.  The  average  explosion 
spectrum  has  a  power  ratio  between  0.5  and  5  Hz  (P^/P^  20  dB 

larger  than  the  corr 3sponding  ratio  for  the  earthquakes,  ex¬ 
cluding  event  229.  bat  this  clear  difference  does  not  occur  for 
the  other  explosion  sites.  The  WESTERN  KAZAKH  presumed  explo¬ 
sions  have  a  spectral  shape  very  similar  to  the  near-by  presumed 
earthquakes  at  Crimea  and  Eastern  Caucacus  (fig. III. 8).  For  the 
NEVADA  explosions  the  large  variation  of  spectra  from  earth¬ 
quakes  in  southwestern  North  America  makes  any  unambiguous  com¬ 
parison  difficult.  Ihe  presumed  explosions  in  WESTERN  RUSSIA 
have  very  strong  high-frequency  energy,  but  no  near-by  earth¬ 
quakes  exist  for  comparison.  The  quakes  at  Mid-Atlantic  ridge 
and  at  western  coast  of  North  America  have  spectra  which  are 
very  deficient  at  high  frequencies  when  compared  with  spectra 
from  events  in  Central  Asia  and  Eastern  Asia  which  are  at  com¬ 
parable  distances,  respectively  (figs  III. 3,  III. 6  and  III. 11). 
The  large  scatcer  and  apparent  increase  of  power  at  high  fre¬ 
quencies  frori  certain  event  regions  (figs  III. 10-12)  is  proba¬ 
bly  due  to  the  very  weak  high-frequency  signal  energy  from  those 
regions.  At  the  high-frequency  end  of  the  observed  spectral 
range  the  assumption  of  noise  stationarity  made  in  the  noise- 
correction  procedure  described  in  app.  2  is  apparently  incorrect 
(gusts  of  wind,  cultural  noise)  and  the  confidence  limits  com¬ 
puted  for  noise  are  too  narrow.  Thus  the  noise  can  accidentally 
influence  estimation  of  weak  high-frequency  signals . 
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111  *8*  Normalized  acceleration  spectra  (asterisks) 
with  averages  (circles).  (1)  The  3  presumed  explosions  in 

3  9  m^LASA^  V??T?RN  KAZAKE  in  the  magnitude  range  3.7  - 
'1  V^ASA)*  p)  A  presumed  explosion  in  the  same  event 
set  with  magnitude  6.0  mb-  The  WESTERN  KAZAKH  events  lie 
an  approximate  distance  of  25  degrees.  (3)  The  2  events 
°f/V"n"  Set  EASTERN  CAUCASUS  with  magnitudes  4.9  and  5  1 

*iSt™Ce  of  30  de^es.  (4)  The  CRIMEA  e™?  with 
magnitude  4.6  and  distance  22  degrees. 


Normalized  acceleration  power  spectra  fa«:+- 
isks)  with  averages  (circles)  m  u  sPectra  (ast 

event  set  GREECE  with  magnitides  3.8  -  4  sT  end 
tance  range  20-28  degrees.  (2)  The  9  auakes\  Jh/  dl: 
set  TURKEY  in  magnitude  range  4  1-5  7  r  the  even' 

range  23-30  decrrepc?  r*\  ‘  ±.  .  mh  anc*  ln  distaj 

nitude  4.9  and^n^nce  18  deaths . “  YUg6s“VIA  wlth  ”» 


Fiq.  III. 10.  Normalized  acceleration  power  spectra  (aster¬ 
isks)  with  averages  (circles).  (1)  The  4  shocks  in  event 
set  ITALY  in  the  magnitude  range  4.8  -  5.2  and  distance 
range  16-18  degrees.  (2)  The  event  in  GIBRALTAR  with  magni¬ 
tude  4.7  and  distance  29  degrees. 
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q  5.0-5. 5  and  in  distance  range  46-61  degrees. 


event  set  NEVADA  In  the  magnitude  range  4.0-S?0*  <2?  ThJ  2 
events  in  the  same  set  with  magnitudes  5.4  and  5  5  This 
event  set  is  at  distance  75  degrees.  <3)3  events  in IS 
Tt^L^t  “TAH/DAKOTA  with  magnitudes  3  8  M  -  4  6  m 

^de^ngeT^I  4inandeat  SV °F  CA™IA’in  ma^' 
above  2-3  Hz 'are 'considered  u^eUaMe^  degr6eS'  Results 


29 


Finally,  some  ratios  between  different  spectra  are  displayed. 
Fig.  III. 13  shows  the  ratio  of  the  average  spectrum  of  6  events 
in  the  set  WEST  COAST  OF  NORTH  AMERICA  to  the  average  spectrum 
of  10  events  in  the  set  CENTRAL  ASIA  (the  anomalous  event  229 
excluded).  Over  the  frequency  range  0.2  to  4  Hz  the  ratio  de¬ 
creases  by  20  dB,  and  the  decrease  has  closely  the  shape  of 
an  exponential  function.  If  the  decrease  is  assumed  to  be  due 
to  a  larger  attenuation  on  the  path  North  America  -  NORSAR 
than  on  the  path  Central  Asia  -  NORSAR,  we  get  for  the  t* 
parameter  a  difference  of  At*  =  0.279  between  the  two  paths  by 
fitting  the  function 

exp  (nfAt*)  (in  ,  i ) 

to  data.  The  corresponding  curve  is  plotted  to  the  figure.  It 
approximate^  the  observations  wich  a  0.5  dB  r.m.s.  scatter. 

The  implied  assumption  is  the  similarity  of  the  two  average 
source  spectra.  All  events  are  within  the  magnitude  range  4.8- 
6.0.  The  spectra  are  smoothed  over  a  bandwidth  of  0.62  Hz. 

Fig.  III. 14  shows  the  ratio  of  the  average  of  6  explosion  spec¬ 
tra  in  the  set  EASTERN  KAZAKH  to  the  average  of  10  earthquake 
spectra  in  the  set  CENTRAL  ASIA.  The  ratio  increases  by  20  dB 
over  the  frequency  range  0.5  to  5  Hz.  A  simple  exponential 
function  which  gave  an  excellent  fit  to  the  previous  ratio  does 
not  fit  now,  even  if  we  disregard  the  ripples  and  consider  only 
the  trend  of  the  ratio.  The  best-fitting  exponential  function 
is  plotted  to  the  figure.  The  best  fitting  function  of  the  form 
ar  where  f  is  frequency  and  a,  m  are  constants  to  be  deter- 
mined  is  also  plotted.  For  m  is  obtained  the  value  1.0.  The  key 
question  is  apparently  whether  fig.  III. 14  means  that  the  source 
spectra  of  earthquakes  and  explosions  are  systematically  differ¬ 
ent  up  to  the  frequency  6  Hz,  i.e.  frequencies  much  higher  than 
the  corner  frequency,  which  for  these  events  with  magnitudes 
above  5.0  is  situated  at  or  below  1.0  Hz. 
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Fig.  III. 13.  Ratio  of  the  average  spectrum  of  earthquakes 
on  west  coast  of  North  America  to  the  average  spectrum  of 
earthquakes  in  Central  Asia.  The  broken  line  describes  the 
best  fitting  exponential  function  and  the  solid  line  the 
best  fitting  power  aM  of  frequency. 


Fig.  III. 14.  Ratio  of  the  average  spectrum  of  Eastern 

e?pl0SiunSl.t0  the  avera9e  spectrum  of  Central  Asian 
earthquakes.  The  broken  lines  describes  the  best  fitting 

exponential  function  and  the  solid  line  the  best  fitting 
power  A  M  of  frequency.  y 
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IV.  EXPERIMENT  WITH  A  DISCRIMINANT 


The  most  striking  variation  among  spectra  from  different  re¬ 
gions  was  found  in  the  relative  levels  of  high  and  low  frequen¬ 
cy  energy.  The  underground  explosion  spectra  also  differ  from 
earthquake  spectra  by  the  relative  deficiency  of  low-frequency 
energy,  when  events  of  equal  magnitude  are  compared.  To  de¬ 
scribe  these  differences  quantitatively,  the  logarithm  of  the 
normalized  third  moment  of  signal  power  about  the  frequency 
ori9in  was  computed  for  all  spectra.  This  quantity  will  be 
called  LTM.  Weichert  used  the  third  moment  of  amplitude  spec¬ 
trum  for  classification,  and  Anglin^"3  added  a  kind  of  noise 
correction  to  his  definition.  The  definition  of  LTM  is 


LTM 


(Z(f)  -  N(f))f3df 


(Z(f)  -  N  ( f )  )  df 


(4.1) 


where  Z(f)  is  an  estimate  of  the  power  spectrum  of  the  time 
interval  containing  the  wanted  signal,  and  N(f)  is  an  estimate 
of  the  background  noise  spectrum.  Thus  the  integrals  in  (4.1) 
contain  no  bias  due  to  noise.  Because  of  the  variances  of  esti¬ 
mates  Z ( f )  and  N(f)  their  difference  can  be  positive  even  when 
the  signal  power  is  nonexistent.  In  app.  2  the  computation  of 
the  confidence  limits  for  deciding  whether  signal  power  is 
really  present  is  described.  The  comparison  of  the  integrated 
powers  to  the  respective  confidence  limits  indicates  whether 
the  LTM  measurement  is  significant. 
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It  turned  out  that  when  the  low  frequency  band  with  largest 
noise  was  rejected  by  defining  the  frequency  limits  as  f^  = 

0.5  Hz  and  f0  =  5  Hz,  the  LTM  values  were  almost  without  ex¬ 
ception  significant  for  the  analyzed  events,  several  of  which 
were  on  the  fringe  of  detection  and  location. 

In  table  VI . 1  20  event  sets  are  listed,  each  from  a  different 

region.  For  each  region  the  average  value  of  LTM  with  standard 
deviation  and  coefficient  of  correlation  with  event  magnitude 
are  shown.  Magnitude  and  distance  information  are  also  given. 

For  the  sets  with  signif icarrtLTM,!^  correlation  parameters 
of  regression  line  fits  to  data  are  given  in  table  VI. 2. 

Comparing  the  geographically  wide  EASTERN  ASIA  set  to  the 
geographically  limited  SAKHALIN,  KURILES  and  SOUTH  OF  HONSHU 
event  sets  it  can  be  observed  that  the  limited  sets  have 
smaller  LTM  standard  deviations  and  larger  LTM,mb  correlation 
coefficients,  and  that  the  average  value  of  LTM  differs  signif¬ 
icantly  among  these  sets.  These  observations  are  explained  by 
regional  variations  of  LTM  values  within  the  wide  EASTERN  ASIA 
set.  These  variations  hide  the  LTM,^  correlation  in  that  set. 

The  regional  averages  of  LTM  are  plotted  to  fig.  IV. 1  as  func¬ 
tion  of  distance.  It  is  apparent  from  the  figure  that  for  earth¬ 
quakes  there  is  no  correlation  with  distance,  that  the  regional 
^■*-fferences  LTM  are  significant  and  that  the  strong  regional 
effects  destroy  the  usability  of  LTM  as  a  classifier,  if  re¬ 
gional  corrections  are  not  included. 

A  major  cause  of  these  regional  variations  is  probably  varia¬ 
tion  of  the  total  anelastic  attenuation  between  the  respective 
ray  paths.  If  the  total  attenuation  parameter  t*  differs  by 
0.1  between  two  ray  paths,  this  will  produce  a  difference  of 
0.2  in  the  LTM  value,  in  average. 
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Table  IV. 2.  Regression  line  fits  to  LTM  values 
as  function  of  magnitude  for  the  regions  with 
significant  correlation  between  the  two  varia¬ 
bles.  N  =  number  of  events,  a  and  b  are  para¬ 
meters  of  the  equation  LTM  =  a  +  bm^.  The  r.m.s. 
scatter  around  the  regression  line  is  also  given, 
(e)  stands  for  presumed  explosions. 


REGION 

N 

a 

b 

r.m.s. 

SAKHALIN 

20 

1.05 

-0.169 

0.109 

KURILES 

12 

2.08 

-0.278 

0.139 

WESTERN  RUSSIA  (e) 

6 

2.15 

-0.210 

0.139 

TURKEY 

9 

0.869 

-0.148 

0.061 

NEVADA  (e) 

6 

1.613 

-0.280 

0.173 
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Fig.  IV. 1.  Regional  averages  of  the  LTM,  a  measure  of  rel¬ 
ative  amount  of  high-frequency  energy,  with  standard  devia¬ 
tions.  Groups  of  presumed  or  known  underground  explosions 
are  plotted  with  filled  circle,  earthquakes  with  open  cir¬ 
cle.  Small  circles  denote  single  earthquakes. 


Tlu>  n  I  (>[)<«  of  the  regression  line  fitted  to  the  LTM,  m.  obsor- 
vaLlonfi  (table  IV. 2)  is  typically  close  to  0.2.  Variation  ol 
magnitude  over  the  whole  observed  range  3.8  -  6.3  wiLl  then 
produce  an  LTM  variation  of  0.5  units.  This  is  only  half  ot 
tie  observed  regional  variation  of  LTM  (of  signal  high-frequen¬ 
cy  contents) .  E.g.  the  average  LTM  values  for  the  event  sets 
MID-ATLANTIC  RIDGE  and  KURILES  are  -0.23  and  0.72  respectively. 

Comparing  earthquakes  with  underground  explosions  in  Eurasia 
it  is  observed  that  explosions  have  generally  higher  LTM  val- 
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QUASI  A 


Fig.  IV. 2.  P  wave  signals  recorded  at  NORSAR  subcirray  06B 
from  a  Central  Asian  earthquake  with  normal  LTM  value 
(quake  A),  from  the  anomalous  event  229  in  Sinkiang  (quake  B) 
and  from  an  explosion  in  Eastern  Kazakh. 


ues  than  earthquakes  at  comparable  distances,  as  expected.  It 
is  difficult  to  decide  whether  this  is  caused  by  regional  va¬ 
riation  of  wave  transmission  capability  at  high  frequencies, 
or  whether  it  is  a  source  effect  because  no  earthquakes  have 
occurred  adjacent  to  explosion  sites.  There  is  a  variation  of 
LTM  values  over  relatively  short  distances,  e.g.  between  Italy 
(low  LTM),  Greece  (high  LTM)  and  Turkey  (low  LTM). 

The  event  229  of  the  CENTRAL  ASIA  event  set  has  an  LTM  value 
of  1.01  which  is  higher  than  the  mean  LTM  for  explosions  in 


Eastern  Kazakh.  Event  229  has  been  located  by  NOAA  to  42.00°N, 
84.58°E,  only  20  km  from  the  location  of  event  108  (41.8°N, 
84.5°E)  which  has  a  perfectly  normal  LTM  value  0.49.  This 
should  exclude  the  possibility  of  anomalously  low  attenuation 
under  event  229.  The  P  wave  signals  recorded  at  a  subarray  of 
NORSAR  from  event  229  are  compared  in  fig.  IV. 2  to  signals  from 
another  Central  Asian  earthquake  and  an  E. Kazakh  explosion.  It 
is  easy  to  see  that  the  signals  of  event  229  resemble  closely 
the  explosion  signals,  and  that  the  high  value  of  LTM  is  really 
describing  the  physical  reality.  The  discriminant  sug¬ 

gests  this  event  to  be  an  earthquake  (at  Nurmijarvi  station  in 
Finland  Mg  =  4.7,  m^NOAA)  =  5.1).  The  existence  of  an  earth¬ 
quake  having  so  anomalous  P  wave  -oectrum  slows  that  in  addi¬ 
tion  to  attenuation  variation,  the  earthquake  spectrum  can  also 
vary  quite  unexpectedly. 

It  is  disturbing  that  the  earthquakes  closest  to  the  explosions 
in  western  and  southwestern  Russia  also  have  high  (explosion¬ 
like)  LTM  values  (event  sets  CRIMEA  and  CAUCASUS) . 

When  the  Nevada  explosions  are  compared  to  earthquakes  in 
western  North  America,  the  astonishinq  fact  appears  that  earth¬ 
quakes  in  the  continental  U.S.  have  relatively  more  high-fre¬ 
quency  energy  (higher  LTM  values)  than  Nevada  explosions.  If 
one  corrects  for  the  effect  of  event  magnitude  by  taking  the 
LTM  to  Increase  by  0.2  units  when  magnitude  decreases  by  1.0 
units,  as  suggested  above,  the  ranges  of  the  LTM  values  for 
the  two  data  sets  become  equal.  The  earthquakes  in  the  Gulf  of 
California  and  at  the  west  coast  of  North  America  behave  well 
giving  lower  LTM  values  than  Nevada  explosions,  in  average. 

Western  North  Timer ica  is  a  tectonically  complex  region.  Lateral 
variation  in  body  wave  transmission  capability  have  been  re¬ 
ported  in  terms  of  P  wave  magnitude  anomalies25'26  and  of  an- 
elastic  body  wave  attenuation  under  seismic  stations2^.  In 
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Fig.  IV. 3.  Observations  on  P  wave  attenuation  in  the  western 
U.S..  Attenuation  above  average  is  deno  -.ed  by  open  symbols, 
attenuation  less  than  average  by  filled  symbols.  Squares  mark 
the  stations  where  attenuation  of  long  period  body  waves  was 
studied^?,  triangles  mark  earthquakes  used  in  P  wave  magni¬ 
tude  anomaly  studies25'26  and  circles  mark  the  earthquakes  for 
which  the  LTM  values  were  computed  in  this  study.  The  asterisk 
denotes  the  Nevada  test  site. 


fig.  IV.  3  those  observations  are  shown  together  with  observa¬ 
tions  made  in  this  study,  plotted  on  a  map  of  the  western  U.S.. 
Observations  indicating  poor  P  wave  transmission  are  shown 
with  open  symbols,  good  transmission  is  indicated  by  filled 
symbols.  A  region  with  poor  P  wave  transmission  runs  in  north- 
south  direction,  comprising  the  Basin  and  Range  region  and 
the  region  around  and  east  of  the  Gulf  of  California.  The 
Nevada  test  site  is  situated  near  the  border  of  this  region. 

As  poor  body  wave  transmission  capability  attenuates  most  the 
high-frequency  waves,  it  is  apparent  why  the  earthquakes  in 
eastern  Utah  and  Dakota  tend  to  have  high  LTM  values  and  the 
Gulf  of  California  shocks  have  low  values.  The  position  of  the 
Nevada  test  site  is  intermediate  to  these,  and  there  is  again 
the  difficulty  of  telling  how  great  actually  is  the  source  ef- 
feet  on  the  LTM  values.  Manchee  observed  spectra  of  Nevada 
and  the  Gulf  of  California  events  and  concluded  that  a  short- 
period  spectral  discriminant  can  be  used  to  separate  explosions 
in  Nevada  from  the  Gulf  of  California  shocks.  Apart  from  any 
improvement  stemming  from  his  particular  definition  of  a  short- 
period  discriminant,  it  seems  that  a  part  of  the  separation 
comes  from  the  strong  attenuation  of  high  frequencies  at  the 
Gulf  of  California  region. 

A  correlation  of  LTM  variation  with  tectonics  appears  when  it 
is  observed  that  high  values  of  LTM  are  associated  with  re¬ 
gions  where  Intermediate  or  deep  earthquakes  occur  in  addition 
fo  shallow  shocks.  Examples  of  this  are  the  event  sets  GREECE 
(average  LTM  0.50),  KURILES  (0.72),  SOUTH  OF  HONSHU  (0.41). 

Low  values  of  LTM  are  associated  with  seismically  active  re¬ 
gions  where  only  shallow  quakes  occur.  Examples  are  TURKEY 
(0.19),  (northern)  ITALY  (0.08),  MID-ATLANTIC  RIDGE  (-0.23), 
WEST  COAST  OF  NORTH  AMERICA  (-0.21),  and  also  SAKHALIN  (0.21). 

This  correlation  can  be  explained  as  caused  by  increased  rigid¬ 
ity  of  the  upper  mantle  in  regions  with  deep  shocks.  The  rigid 
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mantle  can  produce  tensions  necessary  for  earthquakes,  and 
does  not  attenuate  strongly  the  high-frequency  waves.  From  the 
point  of  view  of  plate  tectonics,  one  could  say  that  the  P 
waves  propagate  along  the  cold  and  rigid  plates  through  the 
upper  mantle.  In  case  of  weakly  seismic  regions  or  boundary 
zones  between  seismic  and  aseismic  regions  our  present  know¬ 
ledge  is  not  sufficient  to  predict  the  expected  range  of  LTM 
values.  If  the  upper  mantle  has  a  very  small  attenuation  (high 
rigidity)  under  the  aseismic  regions  where  the  explosions  in 
Eurasia  occur,  the  LTM  (and  TMF )  difference  between  these  ex¬ 
plosions  and  the  earthquakes  of  the  seismic  belt  could  be 
mainly  due  to  a  difference  of  attenuation  in  the  upper  mantle. 

In  figures  IV. 4  (a-c)  LTM  values  of  the  earthquakes  closest  to 
the  explosion  sites  are  plotted  as  a  function  of  distance  from 
the  explosion  site.  The  LTM  values  are  corrected  for  magnitude 
by  using  a  dfLTMJ/dfn^)  of  0.2  and  reducing  the  LTM  values  to 
a  mb  -  earthquake.  The  figures  would  be  essentially  similar 

without  correction.  When  the  earthquake-explosion  site  distance 
decreases,  the  separation  between  earthquake  and  explosion  LTM 
values  also  decreases.  It  seems  that  for  an  earthquake  occurring 
close  to  the  explosion  site,  or  vice  versa,  classification  using 
LTM  measure  would  be  very  uncertain.  Comparison  with  measure¬ 
ments  made  at  Hagfors  station29'30  indicate  that  no  improvement 
is  achieved  by  use  of  the  original  TMF  discriminant. 

It  is  concluded  that  the  LTM  discriminant  is  unsuitable  for 
separating  earthquakes  and  explosions  because  its  separation 
distance  is  small  compared  with  effects  of  regional  attenuation 
variations.  Because  the  TMF  discriminant24  gives  results  which 
are  (after  taking  logarithm  of  TMF)  linearly  correlated  with 
the  LTM  measure  (for  a  sample  of  25  events30  the  correlation 
coefficient  was  0.88  between  the  two  measures)  it  is  expected 
that  the  same  difficulties  are  met  when  using  that  discriminant. 
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respectively  In  fbt  rf-fc+-ar.^  n<*  ri^ht  from  zero  distance 

(Eastern  Kazakh  tist  site^  In  fe?  S?^Uted  fr°m  50,0  N'77-' 
from  )7  i  m  1 1  c  i  K  j  distances  are  computed 

and  earthquake^are^plotted^ith^f illGd"aPdeSUme<^  eXPlosi°- 

spectively.  P  °  with  filled  and  open  symbols, re- 


A  discriminant  employing  short  period  P  wa re  spectrum  should 
either  make  use  of  that  part  of  the  wave  spectrum  which  is  not 
strongly  influenced  by  variations  in  attenuation  (i.e.  fre¬ 
quencies  below  1-2  Hz)  or  should  have  some  kind  of  buiJt-in 
correction  for  attenuation  effects. 


V.  SUMMARY  AND  CONCLUSIONS 

A  large  seismic  array  was  used  for  estimation  of  the  spectra 
of  P  wave  signals  with  proper  corrections  made  for  the  influ¬ 
ence  of  noise  and  instrument  response.  The  spectra  computed  for 
earthquakes  follow  in  principle  the  predictions  from  analytical 
models  of  the  seismic  source.  The  properties  flat  power  spectrum 
at  low  frequencies,  at  hit h  frequencies  an  average  power  de¬ 
crease  proportional  to  oj  j  where  probably  n  =  2  and  between 
these  trends  a  corner  frequency  which  is  associated  with  the 

source  size  agree  well  with  the  spectra  presented  in  chapter  II. 

31 

According  to  Randall  both  the  dislocation  and  stress-relaxa¬ 
tion  models  exhibit  these  features. 

A  five-fold  variation  of  corner  frequency  (0.2  -  1.0  Hz)  for 
shocks  of  magnitude  5.6  is  observed.  If  the  event  229  in  Sinkiang 
with  anomalously  strong  high-frequency  waves  is  interpreted  as 
having  a  very  high  corner  frequency  of  4  Hz,  the  range  of  varia¬ 
tion  is  even  larger  for  events  within  a  rather  narrow  magnitude 

range.  With  an  inverse  proportionality  between  source  size  and 
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corner  frequency  '  this  would  imply  a  strong  variation  in 
earthquake  fault  lenghts . 

The  explosion-generated  P  wave  signal  spectra  agree  in  principle 
with  analytical  models  of  Mueller  and  Murphy,  and  of  Haskell 
(reviewed  and  revised  by  Soggorn  and  UinndfordS  .  The  spectra 
predicted  by  these  models  also  have  a  corner  frequency  and 
an  w  fall-off  at  high  frequencies.  The  main  difference  from 
the  predicted  earthquake  spectrum  is  at  frequencies  below  the 
corner  frequency,  where  the  explosion  power  may  decrease  towards 
zero  frequency,  so  that  a  maximum  appears  at  the  corner  frequen¬ 
cy  .  In  most  cases  the  explosion  P  wave  signals  did  not  have 
measurable  power  at  frequencies  below  the  corner  frequency  due 
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to  noise  interference.  This  was  a  difficulty  in  the  analysis 
of  the  explosion  spectra,  together  with  the  fact  that  the  spec¬ 
tra  at  the  high  corner  frequencies  of  explosions  began  to  be 
distorted  by  the  anelastic  attenuation. 

When  the  corner  frequencies  were  estimated  visually  from  the 
explosion  spectra,  allowing  for  a  slight  effect  of  attenuation, 
they  were  observed  to  lie  at  roughly  three  times  higher  frequen¬ 
cies  than  the  corner  frequencies  of  earthquakes  with  similar 
magnitudes.  This  causes  the  low-frequency  levels  of  the  explo¬ 
sion  spectra  to  remain  below  those  of  the  earthquake  spectra. 

The  difference  may  be  heightened  by  a  decrease  of  explosion  pow¬ 
er  from  corner  frequency  towards  zero  frequency.  At  frequency 
0.3  Hz  the  spectral  power  of  Eastern  Kazakh  explosions  in  mag¬ 
nitude  range  5.2  -  5.8  is  systematically  20  dB  less  than  the 
spectral  power  of  earthquakes  in  Central  Asia  with  similar  mag¬ 
nitudes.  It  is  impossible  to  say  how  this  ratio  develops  at  even 

lower  frequencies  because  the  earth  noisecovers  there  the  explo¬ 
sion  signal. 

Up  to  here  we  have  discussed  a  frequency  range  where  the  effect 
of  anelastic  attenuation  during  wave  transmission  is  of  minor 
importance.  But  in  the  comparison  of  spectra  at  frequencies  above 
2  Hz,  the  attenuation  of  high  frequencies  becomes  a  dominant  fac¬ 
tor.  The  attenuation  variation  within  the  Earth  probably  r*;’ses 
the  systematic  regional  variation  in  the  relative  levels  of  nigh 
and  low  frequency  energy  noted  when  frequencies  up  to  4  or  6  Hz 
are  observed.  This  conclusion  is  supported  by  the  following  ar¬ 
guments:  (a)  In  comparing  the  spectra  from  events  of  same  type 
and  magnitude  but  at  different  regions,  the  shape  of  the  differ¬ 
ence  in  spectra  could  well  be  explained  by  a  different  amount  of 
frequency- independent  Q  type  attenuation  on  the  respective  propa¬ 
gation  paths,  (b)  Spectra  rich  in  high  frequency  energy  were  ob¬ 
served  from  regions  where  in  addition  to  shallow  seismic  activity 
also  deeper  shocks  occur.  Deep  shocks  are  an  indication  of  in- 
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creased  mantle  rigidity  which  is  expected  to  be  coupled  with 
smaller  anelastic  attenuation  of  seismic  waves. 

The  possibility  of  systematic  variation  in  the  earthquake  source 
spectra  between  the  different  regions  should  not  be  neglected, 
though  reasons  for  such  a  variation  are  not  apparent.  The  varia¬ 
tion  of  the  relative  high-frequency  contents  within  source  re¬ 
gions  very  small  in  extent  should  reflect  true  variation  in  the 
source  spectra.  This  was  studied  using  three  event  groups  in 
Eastern  Asia,  each  less  than  a  degree  wide.  The  high-frequency 
contents  of  each  event  was  indicated  by  a  measure  termed  LTM. 
After  an  empirical  magnitude  correction  was  made  to  the  LTM  val¬ 
ues,  their  distribution  around  the  event  group  mean  value  resem¬ 
bled  closely  a  normal  distribution.  But  the  case  is  not  always 
this  simple,  quite  unexpectedly  one  of  the  earthquakes  (229)  in 
Central  Asia  displayed  an  LTM  value  and  spectral  shape  greatly 
from  other  events  in  the  Central  Asian  population.  It 
seems  that  in  addition  to  a  'normal'  source  variation  of  the 
relative  high-frequency  contents  (LTM  value)  within  an  event 
group,  outliers  may  also  appear. 

The  LTM  measure  is  a  logarithm  of  the  noise-corrected  normalized 
third  moment  of  power  in  respect  to  frequency  origin.  It  is  thus 
a  modification  of  the  TMF  discriminant24,  the  motive  for  a  new 
definition  is  to  achieve  a  sound  noise  correction.  Values  of  LTM 
and  log (TMF)  are  linearly  correlated,  so  the  conclusions  made 
for  LTM  apply  also  to  TMF. 

In  comparing  values  of  LTM  for  explosion  and  earthquake  signals 
it  was  observed  that  the  regional  variations  of  both  earthquake 
and  explosion  LTM  values  are  so  large  that  general  comparison 
is  impossible.  It  was  not  possible  to  find  really  comparable 
earthquake  and  explosion  transmission  paths.  It  seems  that  the 
events  should  not  be  separated  by  more  than  1-2  degrees,  for  un¬ 
ambiguous  comparison3^.  The  earthquakes  closest  to  explosion 


sites  also  had  the  most  explosion-like  (high)  LTM  values. 

We  can  roughly  divide  the  present  event  source  regions  to  three 
types:  one  with  only  shallow  seismicity  (low  LTM  values  of 
earthquakes) ,  a  second  with  both  shallow  and  intermediate  or 
deep  seismicity  (high  LTM  values  of  earthquakes)  and  as  a  third 
type  aseismic  platforms  (very  high  LTM  values  for  explosions) . 
In  the  almost  complete  absence  of  absolute  information  on  an- 
elastic  attenuation,  the  possibility  cannot  be  denied  that  the 
LTM  value  differences  also  between  explosions  and  earthquakes 
in  Eurasia  are  mainly  due  to  transmission  path  differences, 
caused  by  a  (hypothetical)  very  low  wave  attenuation  under  a 
stable  platform.  At  least  it  seems  apparent  that  if  earth¬ 
quakes  in  high-attenuation  regions  of  only  shallow  seismicity 
are  compared  to  explosions  on  stable  platforms  for  assessment 
of  the  classification  capability  of  classifiers  like  TMF,  far 
too  positive  results  are  obtained. 

It  is  concluded  that  event  classifiers  using  only  P  waves  can 
best  be  designed  by  using  in  some  suitable  manner  the  one  dif¬ 
ference  between  earthquakes  and  explosions  well  documented  both 
by  theory  and  observation,  the  difference  in  the  respective 
corner  frequencies  of  spectra,  maybe  together  with  the  (less 
well  documented)  behaviour  of  the  spectra  at  frequencies  below 
the  corner  frequency.  It  will  be  necessary  to  take  into  account 
the  predicted  variation  of  corner  frequency  with  event  magni¬ 
tude.  Future  work  will  be  directed  along  these  lines. 
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APPENDIX  1 


SIGNAL -TO-NOISE  RATIOS,  BEAM  LOSSES  AND  STABILITIES 

Estimates  on  P  wave  signal-to-noise  ratio  (SNR)  and  beam  loss 

(due  to  incoherency  of  signals  across  the  array)  as  functions 

of  frequency  have  been  made  for  6  explosions  in  E. Kazakh  and 

9  earthquakes  in  Central  Asia.  The  results  reviewed  here  have 

been  published16.  Events  are  in  the  magnitude  range  5. 0-6.0  mb- 

Power  spectra  of  P  waves  recorded  by  all  single  sensors  of 

NORSAR  were  computed  and  averaged  over  the  sensors,  for  each 

event.  The  average  signal  spectra  were  then  noise-corrected, 

smoothed,  normalized  and  averaged  within  each  group  of  events. 

The  average  earthquake  and  explosion  spectra  were  then  compared 

to  an  average  of  noise  spectra  measured  at  15  random  dates  and 

times  of  year.  Table  1  summarizes  the  SNR  results.  The  column 

marked  f  gives  the  frequency  of  maximum  SNR  and  columns  f,_ 
max  ^  12 

give  frequencies  at  which  SNR  is  12  dB  down  from  the  maximum 


Table  1 .  Summary  of  SNR  measurements 
for  Central  Asia 


event  type 

*12 

f 

max 

f  12 

earthquakes 

0.7  Hz 

2.0  Hz 

4.5  Hz 

explosions 

1.4  Hz 

2.6  Hz 

5.2  Hz 

The  beam  losses  were  also  computed  separately  for  explosions 
and  earthquakes.  For  each  event,  the  signal  power  spectx'um 
of  the  array  beam  was  compared  to  the  average  signal  power 
spectrum  of  the  single  sensors.  Their  ratio  is  the  beam  loss. 
The  losses  were  averaged  within  each  group  of  events.  The  av- 
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cragc  of  the  subarray  beam  spectra  was  similarly  compared  to 
the  average  of  the  single  sensor  signal  spectra  for  each  event, 
and  the  subarray  beam  losses  again  averaged  within  each  event 
group.  The  results  are  shown  in  figs  Al.l  and  A1.2.  The  figures 
show  that  the  signal  losses  approach  the  suppression  of  a  ran¬ 
dom  signal  at  the  frequency  5  Hz  (random  signal  is  expected  to 
be  suppressed  by  21  dB  in  the  array  beam  and  7.8  dB  in  a  sub¬ 
array  beam) .  The  length  of  time  intervals  taken  for  spectral 
analysis  was  6.4  sec. 


relative  i-n  1  °f  Sfgnal  Power  in  array  beamforming 
leipf1  average  signal  power  at  the  single  sensor 

a  functfon  °5  frequency.  The  losses  are  averaged 
er  several  events.  They  are  computed  separately  for  Central 
Asian  earthquakes  (closed  circles  with  standard  deviations) 
and  explosions  (open  circles)  .  aeviauons) 
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FREQUENCY  (CPS) 


relatfvp2^n  Ihl  l0SS  °f  sfgnaJ-  Power  in  subarray  beamforming 
elative  to  the  average  signal  power  at  the  single  sensor 

a,  function  <?f  frequency.  The  losses  are  averaged 
V subarrays  and  several  events.  They  are  computed  sep- 
arately  for  Central  Asian  earthquakes  (closed  circles  with 
standard  deviations)  and  explosions  (open  circle?. 


TIih  oliHoi  V«,1  MNIi  and  beam  loan  wain  u»a.i  |„  ( |„,  t1tj,  llt|l  it, 

stabilities  of  spectra  computed  by  different  forms  of  a.. ay 
processing.  The  following  equations  given  by  Lacoss  and  Kustei 
give  the  stabilities  sp  and  sfa  of  the  single  sensor  average 
signal  spectrum  and  beam  signal  spectrum,  respectively. 


I  5 
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s  =  - - 

4  (r  +  (1  +  £)  ) 


sb  = 


(KT)2 


where  K  is  the  number  of  sensors  used  in  averaging  the  signal 
spectra  or  forming  the  beam,  r  is  twice  the  ratio  of  single 
sensor  signal  power  to  noise  power,  g_1  is  the  beam  loss  and  L 
is  the  number  of  blocks  used  in  estimating  the  noise  power  spec¬ 
trum.  The  stability  of  the  average  of  subarray  beams  is 


N  (~)  2 

4  (f  +  (i +  i» 


where  N  is  the  number  of  subarrays,  M  is  the  number  of  sensors 

in  each  subarray,  e  is  the  average  subarray  beam  loss  and  r 
and  L  are  as  above. 

The  beam  loss  and  SNR  observations  for  Central  Asian  earth¬ 
quakes  were  put  to  the  equations  above,  the  SNR  being  scaled 
suitable  for  a  4.8  magnitude  shock.  The  resulting  stabilities 
as  functions  of  frequency  are  shown  in  fig.  A1.3.  It  is  con¬ 
cluded  that  the  average  of  subarray  beams  has  a  good  stability 
at  a  wide  range  of  frequencies,  processing  does  not  use  an  ex¬ 
cessive  about  of  computing  time  and  the  observed  spectrum  is 
not  distorted  by  more  than  6-7  dB  in  the  worst  case.  It  is 
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FREQUENCY  (CPS) 


Fig.  A1.3.  Stabilities  for  different  forms  of  precessing 
array  data,  computed  from  loss  and  SNR  observations  for  the 
earthquakes  in  Central  Asia,  scaled  to  m.  4.8.  Continuous 
line  is  beamforming,  broken  line  is  spectraforming  and  dotted 
line  is  spectraforming  of  subarray  beams. 


selected  as  the  short  period  array  processing  method  used  in 
this  study. 

The  stabilities  of  the  signal  spectrum  on  long  period  beam  and 
of  the  average  of  LP  single  sensor  spectra  were  also  computed 
using  the  equations  given  above.  It  was  concluded  that  because 
the  LP  beam  loss  is  apparently  small  (not  more  than  a  few  dB) 
beamforming  is  the  better  LP  processing  method  at  all  signal- 
to-noise  ratios. 
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APPENDIX  2 


NOISE-CORRECTION /NORMALIZATION  AND  AVERAGING  OF  SPECTRA 

Z ( f )  is  taken  to  be  the  power  spectrum  of  a  section  of  recording 
which  is  supposed  to  contain  a  transient  signal  and  noise  (f  = 
frequency) .  The  power  spectrum  of  the  pure  signal  is  assumed 
to  be  | S  ( f )  |  .  Y ( f )  is  the  power  spectrum  of  a  section  of  the 
same  recording/  containing  only  noise.  The  noise  is  assumed  to 
be  a  stationary  random  Gaussian  process  with  zero  mean  and  un¬ 
correlated  with  the  signal.  The  noise  power  Y(f)  has  at  each 

frequency,  by  definition,  a  probability  density  distribution 

2 

identical  to  a  multiple  of  the  X  distribution  with  two  degrees 
2  2 

of  freedom  .  The  recordings  from  which  these  sections  are  taken 
may  be  from  a  single  sensor,  a  subarray  beam  or  the  array  beam. 

An  estimate  of  signal  power  is  made  by  subtraction 

P(f)  =  Z  (f )  -  Y (f )  (1) 

The  expected  value  of  P  is^ 

E(P)  =  1  S |  2  (2) 

whore  llu*  Implicit  dependence  on  frequonry  is  not  written,  for 
brev i ty . 

2 

If  | S }  =0  then  E(P)  is  also  zero  and  two  estimates  of  noise 

power  are  being  subtracted  from  each  other.  Because  the  esti¬ 
mates  are  less  than  perfect,  values  of  P  larger  than  zero  can 
occur.  We  want  to  be  on  guard  against  the  situation  where  such 
occurrence  is  taken  to  indicate  signal  power,  while  true  signal 
power  is  zero  or  practically  zero. 
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The  stability  of  the  estimates  of  z  and  Y  can  be  improved  by 
ensemble  averaging  and  frequency  smoothing.  Considering  ensem- 
le  avei aging,  if  there  are  available  N  parallel  channels  re¬ 
cording  the  same  signal  but  different  realizations  of  noise, 
as  is  the  case  with  an  array  station,  the  signal  and  noise 
section  stabilities  can  be  improved  by  averaging  over  these. 

Because  several  noise  sections  at  different  times  can  be  ana¬ 
lyzed  from  each  channel,  the  number  M  of  noise  sections  can  be 
larger  than  N.  Even  though  the  noise  is  here  assumed  to  be 
roughly  an  ergodic  process,  the  ergodicity  is  in  all  probabil- 

ty  so  incomplete  that  one  does  best  by  including  all  channels 
to  noise  averaging  in  equal  proportions. 

The  averages  are 

N  M 

2  =  N  H  zn  (3) 

n~l  m=l 

When  I S |  =  0  the  probability  density  distributions  for  both 

are  multiples  of  X  distributions  with  2N  and  2M  degrees  of 
freedom. 

The  variance  of  P  given  |s|2  =  0,  when  P  is  estimated  as 

Z"Y/  is  the  sum  °f  the  variances  of  Z  and  Y.  These  can  be 
deduced  from  their  X2  density  distributions.  We  get  for  the 
variance  of  p 

,4, 

When  N,M»  1  the  probability  density  distributions  of  the  quan- 
rties  z,  Y  and  P  approach  closely  normal  distributions.  Then 
it  is  easy  to  define  confidence  limits  for  P  given  [si2  =  o 
The  99%  confidence  limit  has  been  used,  or  if 
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(5) 


P  **  2.58  (£  +  fc>*Y 

we  were  99%  confident  in  making  the  right  decision  when  it  was 
concluded  that  signal  power  is  present.  When  P  did  not  fulfil 
inequality  (5),  the  pov.’er  estimate  at  that  frequency  was  re¬ 
jected. 

When  the  array  beam  is  used,  N  =  1  by  necessity.  To  decrease 
the  variance  of  P,  M  can  be  increased  either  by  analyzing  sev¬ 
eral  successive  noise  sections  of  the  beam  recording,  or  by 
analyzing  the  constituents  of  the  array  beam,  the  single  sensor 
or  subarray  beam  recordings,  and  then  similating  the  noise  re¬ 
duction  in  beamforming  by  dividing  their  average  power  by  the 
number  of  sensors  or  subarray  beams  used  in  forming  the  array 
beam,  respectively. 


When  estimating  the  spectra  from  recordings  of  long  period  sen¬ 
sors  presented  in  this  paper,  the  latter  way  was  used.  In  addi¬ 
tion  to  the  signal  section  spectrum  from  the  long  period  array 
beam,  noise  spectra  from  M  single  sensor  recordings  were  com¬ 
puted,  averaged  and  divided  by  N.  Usually  M  was  10,  the  sensors 
of  the  outermost  sensor  ring  of  NORSAR  were  used  to  avoid  any 
local  noise  anomalies. 

When  N  Is  as  small  aB  1  the  probat ility  density  distribut ion  of 
/.  do  v  1  a  I  on  from  the  normal  distribution.  In  computing  (ho  long 
fjoi  Lod  Hpoctra  with  N  a  1 ,  this  deviation  was  neglected  and 
inoq.  (5)  used.  Usually  some  frequency  smoothing  was  made  and 
in  the  summing  of  the  different  frequency  components  the  central 
limit  theorem  was  shaping  the  distribution  of  Z  towards  normal 
distribution. 

Spectra  from  both  long  and  short  period  sensors  were  smoothed 
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over  frequency.  The  purpose  was  twofold:  the  most  rapid  varia¬ 
tions  of  the  spectra  were  not  needed  at  this  level  of  the  ana¬ 
lysis  of  data,  also  smoothing  increases  the  stability  of  the 
noise  estimates  making  possible  to  dig  deeper  to  noise  for  sig¬ 
nal. 

The  variances  of  power  P  given  | S | 2  =  0  at  frequencies  f1#f .  .. 

fk  are  9iven  by  eq. (4 ) .  The  variance  of  the  average  power  PA 
over  K  frequencies  is  then 

°2<v  -ik  +  k'-jjrF- y2(fk)  (6) 

K  k=l 

are  equal,  the  variance  of  the  average  power  is 
decreased  by  a  factor  of  K  relative  to  the  individual  variances. 
Usually  this  is  not  the  case,  however,  the  variance  always  de¬ 
creases  in  smoothing,  with  a  corresponding  narrowing  of  the 
confidence  limits  of  P^.  The  confidence  limits  of  P^  can  easily 
be  estimated  because  of  its  normal  probability  distribution. 

To  average  spectra  from  several,  say  J,  events  each  spectrum 
is  first  normalized.  For  an  event  j  the.  normalizing  factor  is 
computed  as 

CO 

Fj  =  f  Zj(f)  -  Yj (f)  df 
o 

Power  of  noise  ana  signal  spectra  averaged  over  events  is 

r  i  i  <r  zi,£) 

ave  [Zj(f)]  =  i  Z  -i— 

r  i  i  J.  V£> 

ave  [Vjf £)]  -1 - 

j=l  j 

The  average  pure  signal  power  P  is  then 

Pi 
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PA  =  ave(Zj(f))  ~  aveUj(f)) 


with  variance 


where  the  individual  signal  and  noise  spectra  are  assumed  to 
come  from  ensemble  averaging  over  N  and  M  channels  and  to  be 
smoothed  over  K  frequency  components.  If  the  individual  spectra 
resemble  each  other  the  stability  is  increased  in  averaging 
over  events. 


Eq. (4.1)  for  computation  of  the  third  moment  of  power  has  integ¬ 
rals  over  power  in  denominator  and  nominator.  The  variance  a2 
of  the  denominator  integral  (written  as  sum  over  K  frequency 
components)  is 

°2  '  {  +  n  ;  |i  Y  (fk>  2 

The  variance  of  the  integral  in  nominator  is  similar  save  the 
absence  of  factor  f2.  From  the  central  limit  theorem  follows 
that  the  probability  density  distributions  of  the  two  integrals 
arc  normal,  and  confidence  limits  for  both  can  thus  easily  be 
computed  from  the  standard  deviations. 
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APPENDIX  3 


In  the  following  table  the  events  discussed  in  this  report  are 
listed.  They  are  grouped  geographically,  the  groups  containing 
up  to  23  events.  Events  in  Asia  are  listed  first,  then  events 
in  Europe  and  then  events  in  North  America.  Some  larqe  events 
are  separately  added  to  the  end  of  the  list. 


Preceding  page  b'aak 
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